Mon. Not. R. Astron. Soc. 000, 000-000 (0000) Printed 2 February 2008 (MN WF$t style file v2.2) 



Is there Supernova Evidence for Dark Energy 
Metamorphosis ? 



o 
o 



Ujjaini Alam 1,4 , Varun Sahni 1 ' 5 , Tarun Deep Saini 2 ' 6 and A. A. Starobinsky 



3,7 



1 Inter University Centre for Astronomy & Astrophysics, Pune, India 

2 Institute of Astronomy, Madingley Road, Cambridge, UK 

3 Landau Institute for Theoretical Physics, Kosygina 2, Moscow 119334, Russia 
4 ujfaini@iucaa.emet. in 
5 »"« | ™ iQiucaa. ernet. in 



'varum 
6 tarun€ 
7 alstari 



ast.cam.ac.uk 
Handau.ac.ru 



2 February 2008 



> 



o 
6 

c3 



ABSTRACT 

We reconstruct the equation of state w(z) of dark energy (DE) using a recently released 
data set containing 172 type la supernovae without assuming the prior w(z) —1 (in 
contrast to previous studies). We find that dark energy evolves rapidly and metamor- 
phoses from dust-like behaviour at high z (w ~ at z ~ 1) to a strongly negative 
equation of state at present (w < — 1 at z ~ 0). Dark energy metamorphosis appears 
to be a robust phenomenon which manifests for a large variety of SNe data samples 
provided one does not invoke the weak energy prior p + p 0. Invoking this prior 
considerably weakens the rate of growth of w(z). These results demonstrate that dark 
energy with an evolving equation of state provides a compelling alternative to a cos- 
mological constant if data are analysed in a prior-free manner and the weak energy 
condition is not imposed by hand. 
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1 INTRODUCTION 

One of the most tantalizing observational discoveries of the 
past decade has been that the expansion of the universe 
is speeding up rather than slowing down. An accelerating 
universe is strongly suggested by observations of type la 
high redshift supernovae provided these behave as standard 
candles. The case for an accelerating universe is further 
strengthened by the discovery of Cosmic Microwave Back- 
ground (CMB) anisotropics on degree scales (which indicate 
fitotai — 1) combined with a low value for the density in 
clustered matter fl m ~ 1/3 deduced from galaxy redshift 
surveys. All three sets of observations strongly suggest that 
the universe is permeated by a relatively smooth distribu- 
tion of 'dark energy' (DE) which dominates the density of 
the universe (Qde — 2Q, m ~ 2/3) and whose energy momen- 
tum tensor violates the strong energy condition (p + 3p ^ 0) 
so that iode = p/p < —1/3. 

Although a cosmological constant (w = —1) provides a 
plausible answer to the conundrum posed by dark energy, 
it is well known that the unevolving cosmological constant 
faces serious 'fine tuning' problems since the ratio between 
PA = A/8ttG and the radiation density, p r , is already a 
miniscule pa/ pr ~ 10~ 54 at the electroweak scale (T ~ 100 



GeV) and even smaller, pa/ pr ~ 10 123 , at the Planck scale 
(T ~ 10 19 GeV). This issue is further exacerbated by the 
'cosmological constant problem' which arises because the 
A-term generated by quantum effects is enormously large 
Pa > m-pii where mpi ~ 1.2 x 10 19 GeV is the Planck mass 
(Zeldovich 1968; Weinberg 1989). 

Although the cosmological constant problem remains 
unresolved, the issue of fine tuning which plagues A has 
led theorists to explore alternative avenues for DE model 
building in which either DE or its equation of state are 
functions of time. (Following Sahni et al. (2003) we shall 
refer to the former as Quiessence and to the latter as Ki- 
nessence.) Inspired by inflation, the first dark energy mod- 
els were constructed around a minimally coupled scalar field 
(quintessence) whose equation of state was a function of time 
and whose density dropped from a large initial value to the 
small values which are observed today (Peebles & Ratra 
1988; Wetterich 1988). ('Tracker' quintessence models had 
the advantage of allowing the current accelerating epoch to 
be reached from a large family of initial conditions (Cald- 
well, Dave & Steinhardt 1998).) 

Half a decade after SNe-based observations pointed to 
the possibility that we may be living in an accelerating uni- 
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verse, the theoretical landscape concerning dark energy has 
evolved considerably (see the reviews Sahni & Starobinsky 
2000; Carroll 2001; Peebles & Ratra 2002; Sahni 2002; Pad- 
manabhan 2003). In addition to the cosmological constant 
and quintessence, the current paradigm for DE includes the 
following interesting possibilities: 

• Dark energy with w < — 1 (Chiba, Okabe, & Yam- 
aguchi 2000; Caldwell 2002; Mclnnes 2002; Sahni & Shtanov 
2003; Alam & Sahni 2002; Caldwell, Kamionkowski & Wein- 
berg 2003; Carroll, Hoffman & Trodden 2003; Frampton 
2003; Frampton & Takahashi 2003; Singh, Sami & Dadhich 
2003; Johri 2003) 

• The Chaplygin Gas whose equation of state drops 
from w = at high redshifts to w ~ — 1 today (Kamen- 
shchik, Moschella & Pasquier 2001) 

• Braneworld models in which the source for cosmic 
acceleration rests in the gravity sector rather than in the 
matter sector of the theory (Deffayet, Dvali & Gabadadze 
2002; Sahni & Shtanov 2003; Maeda, Mizuno & Torii 2003) 

• Dark energy models with negative potentials 
(Felder et al. 2002; Kallosh et al. 2002; Alam, Sahni & 
Starobinsky 2003) 

• Interacting models of dark matter and dark en- 
ergy (Amendola 2000; Chimento et al. 2003; Hoffman 2003) 

• Modified gravity and scalar-tensor theories 
(Boisseau et al. 2000; Bertolami & Martins 2000; Bartolo 
& Pietroni 2000; Damour, Kogan & Papazoglou 2002) 

• Dark energy driven by quantum effects (Sahni & 
Habib 1998; Parker & Raval 1999) 

• Dark energy with a late-time transition in the 
equation of state (Bassett et al. 2002; Corasaniti et al. 
2003) 

• Unified models of dark energy and inflation (Pee- 
bles & Vilenkin 1999; Copeland, Liddle & Lidsey 2001; 
Sahni, Sami & Souradeep 2002) etc. 

Faced with the current plethora of dark energy scenarios 
the concerned cosmologist is faced with two options: 

(i) She can test every single model against observations, 

(ii) She can take a more flexible approach and deter- 
mine the properties of dark energy in a model independent 
manner. 

In this paper we proceed along route (ii) and demon- 
strate that model independent reconstruction brings us face 
to face with exciting new properties of dark energy. 

Applying the techniques developed in Saini et al. 
(2000); Sahni et al. (2003) to a new data set consisting of 
172 Supernovae from Tonry et al. (2003) and an additional 
22 Supernovae from Barris et al. (2003) we show that the 
DE equation of state which best fits the data evolves from 
w ~ at z ~ 1 to —1.2 <, w < — 1 today. An evolving 
equation of state of DE is favoured by the data over a cos- 
mological constant for a large region in parameter space. 



2 MODEL INDEPENDENT 

RECONSTRUCTION OF DARK ENERGY 

Supernova observations during the previous decade have 
been pioneered by two teams: The High-z Supernova Search 



Team (HZT) (Riess et al. 1998) and the Supernova Cosmol- 
ogy Project (SCP) (Perlmutter et al. 1999). The enormous 
efforts made by these two teams have changed the way cos- 
mologists view their universe. A recent analysis (Tonry et al. 
2003) of 172 type la supernovae by HZT gives the following 
bounds on the cosmic equation of state (at 95% CL) 



-1.48 < w < -0.72, 



(1) 



when the 2dFGRS prior f2 m /i = 0.2 ± .03 is assumed (Per- 
cival et al. 2001). A similar bound 



-1.61 < w < -0.78, 



(2) 



is obtained for a new sample of high-z supernovae by SCP 
(Knop et al. 2003). 1 

These results clearly rule out several DE contenders in- 
cluding a tangled network of cosmic strings (w ~ —1/3) and 
domain walls (w ~ —2/3). However a note of caution must 
be added before we apply (1) or (2) to the wider class of 
DE models discussed in the introduction. Impressive as the 
bounds in (1) & (2) are, they strictly apply only to dark en- 
ergy having a constant equation of state since this prior was 
assumed both in the analysis of the supernova data set as 
well as in the 2dFGRS study (Tonry et al. 2003; Knop et al. 
2003). Aside from the cosmological constant (w = —1), the 
topological defect models alluded to earlier and the sine- 
hyperbolic scalar field potential (Sahni & Starobinsky 2000; 
Urena-Lopez & Matos 2000; Sahni et al. 2003) no viable 
DE models exist with the property w = constant. Indeed, 
most models of dark energy (Quintessence, Chaplygin gas, 
Braneworlds, etc.) can show significant evolution in w(z) 
over sufficiently large look back times. 

In this paper we shall reconstruct the properties of dark 
energy without assuming any priors on the cosmic equation 
of state. (The dangers of imposing priors on w(z) have been 
highlighted in Maor et al. (2002) and several of our subse- 
quent results will lend support to the conclusions reached in 
this paper.) 



2.1 Cosmological reconstruction of w(z) 

Cosmological reconstruction is based on the observation 
that, in a spatially flat universe, the luminosity distance 
and the Hubble parameter are related through the equa- 
tion (Starobinsky 1998; Huterer & Turner 1999; Nakamura 
& Chiba 1999): 



H(z) = 



d L (z) 
dz \ 1 + z 



(3) 



Thus knowing dh we can unambiguously determine the Hub- 
ble parameter as a function of the cosmological redshift. 
Next, the Einstein equations 



H 1 



8nG r 
— IP, 



+ Pde 

4ttG 
3H' 2 



a 4%G \ ^ „ . 

7jp = mt^ + 3p > ) > 



(4) 



1 It is interesting that, when no priors are set on f2 m , the dark 
energy equation of state becomes virtually unbounded from below 
and has a 99% confidence limit of being < — 1 ! (Knop et al. 2003) 
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are used to determine the energy density and pressure of 
dark energy: 

r 2 



PDE — ^critical — Pm — 



8tvG 



(i-n m (x)) 



(5) 



where p C riticai = 3H 2 /8ivG is the critical density of a FRW 
universe. The equation of state of DE w e ff = Pde/pde fol- 
lows immediately (Saini et al. 2000) 



w eff (x) 



2q(x) 



(2x/3) d(lni?)/da;-l 



3(l-n m (a;)) 1 - (H /H)*n 0n 



(6) 



where fiom = 87rGpom/3-ff ( 2 , £ = 1 + 2. In quintessence 
models and in ACDM, the equation (6) determines the true 
'physical' equation of state of dark energy. However the sub- 
script 'eff ' in w e ff stresses the fact that this quantity should 
be interpreted as an 'effective' equation of state in DE mod- 
els in which gravity is non-Einsteinian or in models in which 
dark energy and dark matter interact. Examples of the for- 
mer include Braneworld models and scalar-tensor theories. 
It is well known that in a large class of Braneworld models 
the Hubble parameter does not adhere to the Einsteinian 
prescription (4) since it includes explicit interaction terms 
between dark matter and dark energy (Deffayet, Dvali & 
Gabadadzc 2002; Sahni & Shtanov 2003). In this case the 
equation of state determined using (6) can still be used to 
characterize DE, but physical interpretations of w e g need to 
be treated with caution. 2 

One route towards the meaningful reconstruction of 
w(z) lies in inventing a sufficiently versatile fitting function 
for either cLl(z) or H(z). The parameters of this fitting func- 
tion are determined by matching to Supernova observations 
and w(z) is determined from (3) and (6). 3 Our reconstruc- 
tion exercise will be based upon the following flexible and 
model independent ansatz for the Hubble parameter (Sahni 
et al. 2003) 



H(x) = H [fir, 



,x 3 + A + A lX + A 2 x 2 ] 



(7) 



where x — 1 + z. This ansatz for H(z) is exact for the cos- 
mological constant w = — 1 (Ai = A2 = 0) and for DE 
models with w = -2/3 (A = A 2 = 0) and w = -1/3 
(Ao — Ai = 0). It has also been found to give excellent 
results for DE models in which the equation of state varies 
with time including quintessence, Chaplygin gas, etc. (Sahni 
et al. 2003; Alam et al. 2003). The ansatz (7) is equivalent 
to the following expansion for DE 



PDE = P0c(A + A\X + A 2 x 2 + A 3 x :i ) 



(8) 



2 One way around this difficulty is to define observables solely 
in terms of H and its derivatives (called 'Statcfindcrs' in Sahni 
et al. 2003). A detailed discussion of these issues can be found 
in Alam et al. (2003). 

3 Alternatively one could apply an ansatz to w(z) itself (Chiba & 
Nakamura 2000; Weller & Albrecht 2002; Corasaniti & Copeland 
2003; Gerke & Efstathiou 2002; Maor et al. 2002; Lindcr 2003). 
See Alam et al. (2003) for a summary of different approaches to 
cosmological reconstruction. Non-paramctric approaches are dis- 
cussed in Wang & Lovelace (2001); Huterer & Starkman (2002); 
Saini (2003); see also Daly & Djorgovsky (2003); Nuncs & Lidsey 
(2003). 



where po c = 3Hq /(8irG) is the present day critical density. 
The condition A3 ^ allows pde to mimic the properties 
of dark matter at large redshifts (A3 <C 1 follows from large 
scale structure constraints). From (7) and (8) we find Q, m = 
fiom + ^3, i.e. the value of Q m in (7) can be slightly larger 
than fiom in this case. 

Substituting (7) into the expression for the luminosity 
distance we get 



d L (z) 



c 



l+z 



dx 



A + Aix + A 2 x 2 



(9) 



The parameters Ao,Ai,A 2 are determined by fitting (9) 
to supernova observations using a maximum likelihood 
technique. This ansatz has only three free parameters 
(0, m , Ai, A2) since Ao + Ai+A 2 = 1 — for aflat universe. 
A note of caution: since the ansatz (8) is a truncated Tay- 
lor expansion in x — 1 + z its range of validity is z < few, 
consequently the ansatz-derived H(z) and d^(z) should not 
be used at higher redshifts. 

Note that the weak energy condition for dark energy 
Pde 0, /ode + Pde ^ has the following form for the 
ansatz (7) : 



Ao + Aix + A 2 x > 0, Ai + 2A 2 x > 



(10) 



provided we assume that the Q. m x 3 term in (7) is totally 
due to non-relativistic dark matter and does not include any 
contribution from dark energy. The demand that the WEC 
(10) be satisfied for all x ^ (i.e. in the past as well as in the 
future) requires Ao, Ai, A 2 to be non-negative. However, the 
demand that the WEC (10) be satisfied in the past (x > 1) 
but not necessarily in the future, leads to the somewhat 
weaker constraint 



Ai + 2A 2 > 0, A 2 > 



(11) 



(Models in which pde(z) < for z < and which violate 
the WEC in the future, have been discussed in Felder et al. 

(2002) ; Kallosh et al. (2002); Alam, Sahni & Starobinsky 

(2003) .) 

The presence of the term Q m x 3 in (7) has two important 
consequences: (i) It ensures that the the universe transits to 
a matter dominated regime at early times (z 1), (ii) It 
allows us to incorporate information (available from other 
data sets) regarding the current value of the matter density 
in the universe. This information can be used to perform a 
maximum likelihood analysis with the introduction of suit- 
able priors on fiom- In further analysis we will assume that 
the Q. m x 3 term in (7) does not include any contribution from 
dark energy. 

We have also studied simple extensions of the ansatz (7) 
by adding new terms A-ix^ 1 and A4X 4 . The A-ix -1 term 
allows w(z) to become substantially less than —1, thereby 
providing greater leeway to phantom models. The A4X 4 term 
allows DE to evolve towards equations of state which are 
more stiff than dust (w = 0); its role is therefore complemen- 
tary to that of A-ix^ 1 . Despite the inclusion of these new 
terms, our best fit to the supernova data presented below 
does not change significantly(choosing A-i — 0.0003 and 
A4 = 0.008), which points to the robustness of the ansatz 
(7) for the given data set. 

We should add that our reason for choosing an ansatz 
to fit H(z) rather than some other cosmological quantity 
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Figure 1. The fractional deviation A\og(dL Ho) /log(d,L Ho) be- 
tween actual value and that calculated using the ansatz (7) over 
rcdshift for different models of dark energy with Qom = 0.3. The 
solid lines represent quintessence tracker models for potential V = 
Vo/tfi a , with a = 2 and 4. The dotted lines show the deviation 
for Chaplygin Gas models with k = 1 and 5 (where k is the ratio 
between CDM and Chaplygin gas densities at the commencement 
of the matter dominated epoch). The dot-dashed line represents 
the SUGRA potential, V = (M 4+a /4> a ) exp[± (4>/M Pl ) 2 ], with 
M = 1.6 X 10~ 8 M Pl ,a = 11. The dashed horizontal line repre- 
sents zero deviation from model values, which is true for ACDM, 
and w = — 1/3, w = —2/3 quiessence models. 



was motivated by the fact that the Hubble parameter is 
directly related to a fundamental physical quantity - the 
Ricci tensor, and is therefore likely to remain meaningful 
even when other quantities (such as the equation of state) 
become 'effective'. (This happens for instance, in the case of 
the braneworld models of dark energy discussed in Deffayet, 
Dvali & Gabadadze (2002); Sahni & Shtanov (2003).) 

The rationale for choosing a three parameter ansatz for 
H (z) is the following. The observed luminosity distance de- 
termined using type la supernovae is rather noisy, therefore 
in order to determine the Hubble parameter from Dl(z) 
and following that the equation of state, one must take 
two derivatives of a noisy quantity. This difficulty can be 
tackled in two possible ways: (i) either one smoothes the 
data over some interval At (binning is one possibility), or 
(ii) we may choose to smooth 'implicitly' by parameterizing 
H(z) through an appropriate fitting function. The number 
of free parameters JV in the fit to H(z) will be related to the 
smoothing interval Az through Az — z max /A/\ Increasing 
N implies decreasing Az which results in a rapid growth of 
errors through AH(z) oc (Az) _3/2 , and Aw(z) oc (Az)~ 5/2 
(Tegmark 2002), therefore in order not to loose too much 
accuracy in our reconstruction we considered 3 parameter 
fits for H(z) in our paper (these correspond to 2 parameter 
fits for w(z)). 

We now test the usefulness of the ansatz (7) in re- 
constructing different dark energy models. The ansatz re- 
turns exact values for ACDM, and w — —1/3, w = —2/3 
quiessence models. In figure 1 we show the accuracy of 
the ansatz (7) when applied to several other dark energy 
models such as tracker quintessence, the Chaplygin gas and 



super-gravity (SUGRA) models. We plot the deviation of 
log(di,i/o) (which is the measured quantity for SNe) ob- 
tained with the ansatz (7) from the actual model values. 
Clearly the ansatz performs very well over a significant red- 
shift range for Qom = 0.3 (Also see appendix B). In fact, in 
the redshift range where SNe data is available, the ansatz 
recovers these models of dark energy with less than 0.5% 
errors. However it would be appropriate to add a note of 
caution at this point. Although figure 1 clearly demonstrates 
the usefulness of the ansatz for some DE models, its perfor- 
mance vis-a-vis other models of DE is by no means guaran- 
teed. By its very construction the ansatz (7) is expected to 
have limitations when describing models with a fast phase 
transition (Bassett et al. 2002) as well as rapidly oscillating 
quintessence models (Sahni & Wang 2000). (The ansatz (7) 
can give reasonable results even for these models provided 
the resulting DE behaviour is suitably smoothed.) For this 
reason, although the bulk of our analysis will be carried out 
using (7) , we shall supplement it when necessary with other 
fitting functions, which will provide us with an independent 
means with which to test the robustness of our reconstruc- 
tion exercise. 

Methodology : 

For our primary reconstruction, we use a subset of 172 
type la Supernovae, obtained by imposing constraints Ay < 
0.5 and z > 0.01 on the 230 SNe sample, as in the primary 
fit of Tonry et al. (2003). For the ansatz (9), we require 
to fit four parameters: (Ho, £lom, Ai, A 2 ). We may use prior 
information on Ho (Ho = 72 ± 8 km s _1 Mpc -1 , Freedman 
et al. (2001)) and fi m (Q mh = 0.2 ± 0.03, Percival et al. 
(2001)). 4 

The measured quantity for this data is y —< 
\og(d,LHo) >, therefore the likelihood function is given by 



C 



A/"exp 



X_ 

2 



(12) 



^2 _ f Vi - yflt(go,Oom,^i,^2) \ 2 ^ 

1=1 \ a% ' 

where M is a normalisation constant. Therefore, the 
probability distribution function in the four-space 
(H , SJom, A-i, A 2 ) is 



P(H , Q m,A 1 ,A 2 ) oc exp 



X_ 

2 



Pr(n m/i)Pr(ffo) • (14) 



where Pr refers to the priors applied on the parameters of 
the system. 

Our goal is to reconstruct cosmological parameters such 
as the equation of state w(z) — w(z; Qom, Ai, A2), therefore 
we marginalise over Ho and obtain the probability distribu- 
tion function in the (fiom, A\, A2) space: 



P(Q m,A 1 ,A 2 )= / P(HoM 0m ,A 1 ,A 2 )dHo. 



(15) 



4 One should note however that the prior on Qom is not model 
independent since it relies on the ACDM model to project from 
rcdshift space to real space. Results coming from the use of this 
prior should therefore not be taken too literally in the present 
context. See Kunz et al. (2003) for an interesting discussion of 
related issues. 
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Figure 2. The (A\, A2) parameter space for the ansatz (7). The 
light grey shaded area shows the allowed region if dark energy sat- 
isfies the weak energy condition both currently and in the past: 
w(z) ^ —1,2 ^ 0. The x 2 surface has two minima, a shallow 
minimum at Aj = 0.177, A 2 = -0.119 with xLallow = 1-040 2 
and a deeper minimum at Ai = —4.360, A2 = 1.829 with Xd eep = 
1.0056. The deeper minimum is marked by a bullet. The solid con- 
tours surrounding the deeper minimum are la, la, 3a contours of 
constant A\ 2 where Ax 2 = X 2 — X 2 ecp 1 Similarly the dashed con- 
tours surrounding the shallower minimum are la, 2a, 3<r contours 
of constant Ax 2 where Ax 2 = X 2- Xghaiiow- Note that the ACDM 
model (marked by a solid star) corresponds to A\ = A2 = which 
is very close to the shallow minimum. 



In order to do this, we have to define the bounds of a four- 
dimensional volume in (Ho, fiom, Ai, A2). The bounds of Ho 
are taken at 5a of the HST prior. For fiom, t ne natural 
choice is ^ Qom ^ 1. It is not immediately obvious what 
the bounds should be for Ai,A2. We choose a sufficiently 
large rectangular grid for Ai , A2 (roughly corresponding to 
— 6 <, w < 5) which includes most known models of dark en- 
ergy. This bound is merely a matter of convenience and does 
not affect our results in any way. After marginalisation, we 
have a grid in (fiom, Ai, A 2 ) space on which P(f2 m, Ai, A 2 ) 
is specified at each point. We may now proceed in two ways. 
Firstly, we may choose to fix fiom at a suitable constant value 
(e.g., Q,om = 0.3) thereby obtaining a grid in the (Ai,A 2 ) 
plane with P (the probability if f2om is known to be an ex- 
act value) denned at each point. For a particular redshift, 
we may then calculate w(z; fiom, Ai, A2) at each point of 
the grid. This would yield results that would hold true if 
Oom were known exactly. Instead of using the exact value 
of fiom, we may use the prior information about it available 
to us (Slomh — 0.2 ±0.03), and calculate w(z; Hom, Ai, A2) 
at each point of a three-dimensional grid, the probability 
P at each point being known. Therefore, at any given red- 
shift z, u>(fiom, Ai, A2) can be tagged with a numerical value 
P(f2om, Ai,A2). Starting from the best-fit w(z) (the value 
at the peak of the probability distribution), we may move 
down on either side till 34% of the total area is enclosed 
under the curve, thus obtaining asymmetric la bounds on 
w(z). The 2a, 3a bounds can be similarly obtained. 
Results : 
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Figure 3. The deviation of H 2 /H 2 from corresponding ACDM 
values over redshift for the ansatz (7). The thick solid line shows 
the best-fit, the light grey contour represents the la confidence 
level, and the dark grey contour represents the 2a confidence level 
around the best-fit. The dashed horizontal line denotes ACDM. 
fi()m = 0.3 is assumed. 



We first show preliminary results for which the matter 
density is fixed at a constant value of Qom = 0.3. A detailed 
look at the \ 2 surface in the (A\, A2) plane (figure 2) reveals 
the existence of two minima in x 2 , a shallower one close 
to ACDM (Ai = 0.177, A 2 = -0.119, w = -1.03, x 2 = 
1.0402), and a deeper minimum at Ai — —4.360,^42 = 
1.829, w = -1.33, x 2 = T0056. Wc would like to draw 
the readers attention to the fact that imposing the prior 
w(z) ^ — 1 (z ^ 0) amounts to disallowing a significant re- 
gion of parameter space (the unshaded region in figure 2). 
Consequently an analysis which assumes w(z) —1 loses 
all information about the region 2a around the deeper mini- 
mum ! Since we have no reason (observational or theoretical) 
to favour either minimum over the other, we shall always 
choose the deeper minimum as our best-fit in all the subse- 
quent calculations. 

In the figure 3, we plot the deviation of the squared 
Hubble parameter H 2 /Hq from ACDM over redshift for the 
best-fit. We note that the quantity H 2 /H 2 has a simple 
linear relationship with the parameters of the fit (Eq 7), 
therefore the errors in this quantity increase with redshift. 
Another quantity of interest is the energy density of dark 
energy. For this ansatz, p D E = Pde/poc = A + Aix + A 2 x 2 
(where poc = 3H / 8irG is the present day critical density) . 
The figure 4 shows the logarithmic variation of pde with red- 
shift. In this figure too the errors increase with redshift. An 
interesting point to note is that initially, dark energy den- 
sity decreases with redshift, showing the phantom-like na- 
ture (w < — 1) of dark energy at lower redshifts of z < 0.25, 
while at higher redshifts, the dark energy density begins to 
track the matter density. Before moving on to the second 
derivative of the luminosity distance (e.g., the equation of 
state) we may obtain more information from the dark energy 
density by considering a weighted average of the equation of 
state : 

1 + ^Aln(TTI)/ (1 + W(2)) T$i' ( 16 > 
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Figure 4. The logarithmic variation of dark energy density 
Pde/P0c (where poc = 3H^/8ttG is the present day critical den- 
sity) with redshift for the ansatz (7). The thick solid line shows 
the best-fit, the light grey contour represents the la confidence 
level, and the dark grey contour represents the 2a confidence 
level around the best-fit. The dashed horizontal line denotes 
ACDM and the dotted line represents matter density f2o m (l+z) 3 , 
f^Om = 0.3 is assumed. 



Table 1. The weighted average w (eq 16) over specified redshift 
ranges. The best-fit value and la and 2a deviations from the 
best-fit arc shown. 



la 



2a 



- 0.414 



-0.969 



+0.120 +0.198 
-0.089 -0.199 



0.414-1 -0.108 i-° 2 360 

1 - 1.756 0.069 ±°;J°° ±S:ifS 



where A denotes the total change of the variable between 
integration limits. This quantity can be elegantly expressed 
in terms of the difference in energy densities over a range of 
redshift as 



1 + w = 



1 A InpDE 
3 A ln(l + z) 



(17) 



Thus the variation in the dark energy density depicted in 
figure 4 is very simply related to the weighted average equa- 
tion of state ! 

Table 2. x P er degree of freedom for best-fit and ACDM models. 
wo is the present value of the equation of state of dark energy in 
best-fit models. Xp r refers to the best fit after imposition of the 
WEC prior w{z) -1 (z ^ 0). 



Best-fit 



Confidence levels w(z) > -1 ACDM 



f^Om 


wo 


X 2 ■ 

A min 


Xl CT 


X2ct 


Xp r 


x 2 


0.10 


-1.093 


1.0077 


1.0213 


1.0442 


1.0359 


1.1242 


0.20 


-1.198 


1.0071 


1.0207 


1.0436 


1.0384 


1.0663 


0.30 


-1.334 


1.0056 


1.0192 


1.0421 


1.0409 


1.0417 


0.40 


-1.470 


1.0043 


1.0179 


1.0408 


1.0578 


1.0638 


0.50 


-1.606 


1.0038 


1.0174 


1.0403 


1.0912 


1.1168 



In table 1 we show the values of to obtained using differ- 
ent ranges in redshift for our best-fit with corresponding la 
and 2a errors. We have taken the ranges of integration to be 
approximately equally spaced in ln(l + z), with the upper 
limit set by the furthest supernova known at present. The 
values of w may be calculated using the equation (16) (which 
uses the second derivative of the luminosity distance), or 
they can simply be read off from figure 4 using equation (17) . 
From this table, a "metamorphosis" in the properties of dark 
energy occurring somewhere between z ~ and z ~ 1 can be 
clearly seen (note that, effectively, one needs to differentiate 
<1l{z) only once to come to this conclusion). 

We now reconstruct the equation of state of dark energy 
which, for the ansatz (7), has the form 



w(z) = -1 + 



Ai + 2A 2 x 



X = 1 + z 



(18) 



3 A + A x x + A 2 x 2 ' 

(Note that, since w(z) was derived from the ansatz (7), its 
domain of validity is z < few.) In figures 5 (a), (b), (c), 
we show the evolution of the equation of state w(z) with 
redshift for different values of Qom- The la and 2a limits 
are shown in each case. The % 2 P er degree of freedom for 
the best-fit for the different cases is given in Table 2. We 
find that for 0.2 < Slom < 0.4, the behaviour of the equation 
of state does not change significantly with change in the 
matter density. However, for larger values of Oo m , a smaller 
current value of w = w(z — 0) is preferred. In all three 
cases considered, the present value of the equation of state 
is w o <, — 1-2 for the best-fit, and the equation of state rises 
steeply from w < — 1.2 to w ~ with redshift. In fact, the 
behaviour of w appears to be extremely different from that 
in ACDM (w — —1). We note here that, for this analysis, 
the errors on w appear to decrease with redshift. This may 
appear counter-intuitive, since there are fewer SNe at higher 
redshifts, but this is merely a construct of the fact that w 
depends non-linearly on the parameters of the ansatz (see 
appendix A). 

Quintessence models satisfy the weak energy condition 
(WEC) p+p = <j> 2 > and it would be interesting to see how 
the imposition of the WEC as a prior on the equation of state 
will affect the results of our analysis. We therefore perform 
the same analysis as above with the added constraint wo ^ 

— 1 (note that this implies w(z) > —1 for all z > for 
our fitting function of H(z) provided A 2 > 0). The results 
are shown in the figures 6 (a), (b), (c). We see that in this 
case the errors are larger and the evolution of the equation 
of state with redshift follows a much gentler slope. Such 
an equation of state would be largely consistent with the 
cosmological constant model. (These results are in broad 
agreement with an earlier analysis of Saini et al. (2000) in 
which a smaller SNe data set was used and a different ansatz 
for the luminosity distance was applied.) 

In Table 2, we show how the x 2 for the best-fit evolving 
dark energy models compare with that for ACDM. We find 
that Xacdm > Xbcst-flt always. For D. 0ln = 0.3, the value of 
Xacdm is just within 2a of the best fit x 2 , but for fi 0m = 
0.2, or fiom = 0.4, Xacdm is outside the 2a limits of the 
best-fit. It is also noteworthy that when the prior w(z) > 

— 1 (z ^ 0) is used, the best-fit model has a slowly evolving 
equation of state with wo = — 1 and the x 2 for the best-fit 
becomes smaller for a smaller value of the matter density. 
When no priors are assumed on w, the trend reverses, and 
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(a) (b) (c) 

Figure 5. The evolution of w(z) with redshift for different values of Qom- The reconstruction is done using the polynomial fit to dark 
energy, equation (7). In each panel, the thick solid line shows the best-fit, the light grey contour represents the Icr confidence level, and 
the dark grey contour represents the 2cr confidence level around the best-fit. The dashed line represents ACDM. No priors are assumed 
on w(z). The x 2 per degree of freedom for each case is given in Table 2. 
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Figure 6. The evolution of w(z) with redshift for different values of Qom, using the prior w(z) — 1,2 ^ 0. The reconstruction is 
done using the polynomial fit to dark energy, equation (7). In each panel, the thick solid line shows the best-fit, the light grey contour 
represents the Icr confidence level, and the dark grey contour represents the 2a confidence level around the best-fit. The hatched region 
is forbidden by the prior w(z) ^ —1. The dashed line represents ACDM. The \ 2 per degree of freedom for each case is given in Table 2. 



better fits are obtained for larger values of f2om- From this it 
appears that at least at la the evolving dark energy model 
is favoured over ACDM, and it does as well, if not better at 
the 2<j level, depending upon the value of the present-day 
matter density. 

Using Priors on fio m : 

Instead of assuming an exact value for fiom, which is 
somewhat optimistic given the present observational sce- 
nario, we may use the 2dF prior on Q 0rn and calculate w(z) 
as a function of (fiom, Ai, A2). It should be noted here that 
the 2dF error bars on fiomfo have been calculated using 2dF 
data in conjunction with CMB, and this calculation assumes 
a ACDM model, therefore this prior should be used more as 
a benchmark for the value of fiom rather than as an absolute 



when considering evolving dark energy models. The resul- 
tant "marginalised" w is shown as a function of the redshift 
in figure 7 (a). The nature of the equation of state for the 
analysis with the added prior w(z) ^ — 1 (z ^ 0) is shown in 
figure 7 (b). We find that the general nature of evolution of 
the equation of state is not changed by adding this extra in- 
formation on the matter density. If no priors are assumed on 
the equation of state to begin with, w(z) still rises sharply 
from w < — 1 up to to ~ at maximum redshift and the 
analysis appears to favour a fast-evolving equation of state 
of dark energy over the standard ACDM model. If a prior 
to ^ —1 is assumed, then the marginalised equation of state 
is more consistent with the cosmological constant. From this 
we see that marginalisation over flom does not lead to any 
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Figure 7. Evolution of w(z) with redshift using the 2dF prior Qomh = 0.2 ± 0.03, with (a) no priors on w(z), and (b) the prior 
w(z) —1, z ^ 0. The reconstruction is done using the polynomial fit to dark energy, equation (7). In both panels, the thick solid 
line represents the best-fit, the dashed line represents ACDM, the light grey contour represents la confidence level, and the dark grey 
contour represents the 2a confidence level. In the right hand panel, the hatched region is forbidden by the prior w(z) — f. 



significant change in our results. In the subsequent sections, 
we will show our results for Qom =0.3. 

From the above analysis, we find, therefore, that our 
results change significantly depending upon whether or not 
the prior w > —1 is imposed. We saw earlier that in the 
absence of any prior on w(z), the best-fit equation of state 
rose from w < — 1 at z = to w ~ at z ~ 1. By imposing 
a prior on the equation of state, we effectively screen off 
a sizeable part of the parameter space (see figure 2), and 
therefore the reconstruction is forced to choose its best-fit 
away from the true minima of the \ 2 surface. The effect of 
imposing a prior on w(z) is therefore to make the best-fit 
w(z) grow much more slowly with w = —1 being preferred 
at 2 = 0. Our results show that the reconstructed equation 
of state with the prior w ^ — 1 is in good agreement with a 
cosmological constant at the 68% CL. However, if no prior is 
imposed, then the steeply evolving dark energy models are 
favoured over the cosmological constant at la, and are at 
least as likely as the cosmological constant at the 2a level. 

Age and Deceleration Parameter of the Uni- 
verse: 

We may also use this ansatz to calculate other quantities 
of interest, such as the age of the universe, t(z), and the 
deceleration parameter, q(z): 



t(z) 
«(*) 



dx 



l + z 

a 



xH(x) 



E _i 



(19) 
(20) 



where x = 1 + z. 

In figure 8 we plot the evolution of the age of the uni- 
verse with redshift. We find that the best-fit age of the uni- 
verse today is to = 12.8 Gyrs if the Hubble parameter is 
taken to be Ho = 72 km s _1 Mpc -1 , which is slightly lower 
than the age of a ACDM universe, to = 13.4 Gyrs (both 
values are for Q m = 0.3). At the 2a level, the age of the 
universe today would vary between 11.2 < to ^ 13.6 Gyrs. 




Figure 8. The age of the universe, Hot(z), is shown as a function 
of the redshift for Qom = 0.3. The reconstruction is done using 
the polynomial fit to dark energy, equation (7). The thick solid 
line represents the best-fit, the light grey contour represents la 
confidence level, and the dark grey contour represents the 2a 
confidence level. The upper (lower) dashed line represents ACDM 
(SCDM). 



Figure 10 shows the evolution of the deceleration pa- 
rameter with redshift. We find that the behaviour of the 
deceleration parameter for the best-fit universe is quite dif- 
ferent from that in ACDM cosmology. Thus, the current 
value of qo — —0.9 is significantly lower than qo — —0.55 
for ACDM (assuming fiom = 0.3). Furthermore the rise of 
q(z) with redshift is much steeper in the case of the best-fit 
model, with the result that the universe begins to acceler- 
ate at a comparatively lower redshift z ~ 0.3 (compared 
with z — 0.7 for ACDM) and the matter dominated regime 
(q ~ 1/2) is reached by z ~ 1. 

Using Priors on Age of the Universe : 
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H = 72 ± 8 km s" 1 Mpc" 1 




(d) (e) (f) 

Figure 9. The evolution of w(z) with redshift for Qom = 0.3 using a Gaussian prior on the age of the universe today: to = 13 ± lGyrs 
(figures (a) and (d)), to = 14 ± lGyrs (figures (b) and (e)), and to = 15 ± lGyrs (figures (c) and (f)). The reconstruction is done using 
the polynomial fit to dark energy, equation (7). The top panel is obtained using a Gaussian prior Ho = 72 ± 8 km s" 1 Mpc -1 , while for 
the bottom panel, Ho = 66±6 km s — 1 Mpc -1 . No priors arc assumed on w(z). In each panel, the thick solid line shows the best-fit, the 
light grey contour represents the Icr confidence level, and the dark grey contour represents the 2a confidence level around the best-fit. 
The dashed line represents ACDM. The \ 2 P er degree of freedom for each case is given in Table 3. 



Table 3. \ 2 P er degree of freedom for the ansatz (7) which best 
fits the SNe data after different age priors are imposed, wq is the 
present value of the equation of state of dark energy in best-fit 
models. 



H 


to 


Best-fit 


km s -1 Mpc -1 


Gyrs 


wo 


X 2 ■ 

^min 


72 ± 8 


13 ± 1 


-1.271 


1.0062 




14 ±1 


-1.099 


1.0197 




15 ± 1 


-0.904 


1.0407 


66 ± 6 


13 ± 1 


-1.553 


1.0139 




14 ±1 


-1.324 


1.0057 




15 ± 1 


-1.153 


1.0146 



Important consistency checks on our best-fit Universe 
may be provided by observations of the age of the Uni- 
verse. Unfortunately, estimates of the age of the universe 
from different methods can produce widely varying results 



one reason for which is that estimates of the Hubble param- 
eter itself can vary significantly. For instance, the HST key 
project yields Ho = 72 ± 8 km s -1 Mpc -1 , while studies of 
the Sunyaev-Zeldovich effect in galaxy clusters give a sig- 
nificantly lower value Ho = 60 ± 10 km s -1 Mpc -1 (Krauss 
2001). Estimates of the ages of the oldest globular clusters 
suggest to = 12.5 ± 2.5 Gyrs, at the 95% confidence level 
(Krauss 2001; Krauss & Chaboyer 2001; Gnedin, Lahav & 
Rees 2001; Hansen et al. 2002; Gratton et al. 2003; Marchi 
et al. 2003) and this age estimate is consistent with sev- 
eral other measurements including observations of eclipsing 
spectroscopic binaries (Thompson et al. 2001; Chaboyer & 
Krauss 2002), results from radioactive dating of a metal- 
poor star (Cayrel et al. 2001) and WMAP data (Spergel 
et al. 2003) (see also Alcaniz, Jain and Dev (2002)). The 
results from the WMAP experiment suggest to = 13.4 ± 0.3 
Gyrs with a Hubble parameter Ho — 72 ±5 km s -1 Mpc" 1 , 
for ACDM cosmology (which satisfies the WEC). Adding 
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Figure 10. Evolution of deceleration parameter of the universe, 
q(z) with rcdshift for Qom = 0.3. The reconstruction is done 
using the polynomial fit to dark energy, equation (7). The thick 
solid line represents the best-fit, the light grey contour represents 
la confidence level, and the dark grey contour represents the 2<r 
confidence level. The dashed line represents ACDM. 



SDSS and SNe la data to WMAP, Tegmark et al. (2003) 
find an age of t = U.lH Gyrs for a slightly closed ACDM 
universe with Ho = 66|;I km s _1 Mpc -1 . Although these 
results cannot be carried over to evolving dark energy mod- 
els including those implied by our best-fit reconstruction 
(which violate the WEC) they provide an indication of the 
range within which the age of the universe might vary. Keep- 
ing in mind these various results, we use two different pri- 
ors on the Hubble parameter: Ho — 72 ± 8 km s _1 Mpc -1 
(1(7 bound from HST; Freedman et al. (2001)), and H = 
66 ± 6 km s _1 Mpc -1 (approximate bound from WMAP, 
SDSS, SNe la; Tegmark et al. (2003)). For each case, we 
choose three different Gaussian priors on the present age of 
the universe: to = 13 ± 1 Gyrs, 14±1 Gyrs, and 15 ± 1 Gyrs 
respectively, and perform the reconstruction for a Qom = 0.3 
universe. The results are shown in the figure 9. We find that, 
for a Hubble parameter of H = 72 ± 8 km s _1 Mpc -1 , and 
with an additional prior on the age of the universe to = 13±1 
Gyrs, the best-fit remains nearly the same, showing a rapid 
evolution of the equation of state from w ~ at z ~ 1 to 
w ~ —1.2 at z = 0, and the errors become narrower. As the 
age is increased, the best-fit equation of state evolves more 
slowly, and the Xdof a l so increases (see Table 3). For the 
prior Ho = 66 ± 6 km s _1 Mpc -1 , wc find that the lowest 
Xdof is obtained for the age prior of to = 14 ± 1 Gyrs, which 
once again matches our best-fit. It should be noted that the 
errors are smaller in all cases, even though the \ 2 m ay be 
larger. We must remember that the addition of a new prior 
which is consistent with the underlying dataset would lead 
to a natural reduction in errors. However, the addition of a 
prior inconsistent with the dataset would lead to a shift of 
the likelihood maximum as well as a reduction in errors, and 
the results would then fail to reflect the actual information 
present in the dataset. That this is happening here for the 
higher values of age can be seen from the fact that although 
the errors are reduced, the Xdof is actually larger. Therefore 
priors from other observations should be added prudently to 



ensure that they do not lead to incorrect representation of 
the data. Since there is as yet no clear model independent 
consensus on the age of the universe, the results we obtain in 
this section should be interpreted with a degree of caution. 

Figure 10 shows the evolution of the deceleration pa- 
rameter with redshift. We find that the behaviour of the 
deceleration parameter for the best-fit universe is quite dif- 
ferent from that in ACDM cosmology. Thus, the current 
value of qo — —0.9 is significantly lower than go — —0.55 
for ACDM (assuming flom = 0.3). Furthermore the rise of 
q(z) with redshift is much steeper in the case of the best-fit 
model, with the result that the universe begins to acceler- 
ate at a comparatively lower redshift z — 0.3 (compared 
with z ~ 0.7 for ACDM) and the matter dominated regime 
(q ~ 1/2) is reached by z ~ 1. 



2.2 Robustness of Results : 

Based on the above analysis, it is tempting to conclude 
that the dominant component of the universe today is dark 
energy with a steeply evolving equation of state which 
marginally violates the weak energy condition. (Of course, 
the less radical possibility of weakly time dependent dark 
energy satisfying the weak energy condition remains an al- 
ternative, too.) However, before any such dramatic claims 
are made, we need to check if our results are in any fashion 
a consequence of inherent bias in the statistical analysis it- 
self, or in the sampling of the data. We therefore perform 
the following simple exercises to satisfy ourselves of the ro- 
bustness of our results. 

Using Different Subsets of Supernova Data : 
In an attempt to understand how the nature of the re- 
constructed equation of state is dependent on the distri- 
bution of data, we perform the reconstruction exercise on 
different samples of data. We have confined ourselves to 
the case where Oo m = 0.3 for these exercises. Firstly, wc 
may exclude the SCP data points from the 172 SNe pri- 
mary fit, leading to a subsample of 130 SNe. We call this 
the HZT sample. Figures 11(a) and 12(a) show the result 
of performing the analysis on this subsample without any 
constraints. The \ 2 P er degree of freedom for the best-fit 
is Xhzt — 0.9707, which is lower than Xacdm — 0.9939 for 
this sample. In this case we find that, though the error bars 
are slightly larger, overall the dark energy density behaves 
in the same way as before (compare figure 11(a) with fig- 
ure 4), showing phantom like (w < —1) behaviour at lower 
redshifts and tracking matter at higher redshifts. The equa- 
tion of state of dark energy also evolves much in the same 
way as when the entire sample is used (compare figure 12(a) 
with figure 5(b)), starting at uio < — 1 and evolving rapidly 
to w ~ 0. We may also use a sample complementary to this 
sample, where all the SCP data points published till date are 
considered, along with the low redshift Calan-Tololo sample. 
This leads to a sample of 58 SNe (Perlmutter et al. 1999; 
Knop et al. 2003), which we call the SCP sample. Using this 
sample, we obtain the figures 11(b) and 12(b). The best-fit 
has a chi-squared per degree of freedom: Xscp ~ 1.0147, 
lower than Xacdm — 1-0369 for this sample. We find that 
here too, the dark energy density initially decreases and then 
starts tracking matter. The equation of state shows signs of 
rising steeply at low redshifts, but since the highest redshift 
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Figure 11. The logarithmic variation of dark energy density Pde( 2 )/P0c (where poc = 3Hq/8wG is the present day critical density) 
with redshift for Uom = 0.3, using (a) HZT data, and (b) SCP data. The reconstruction is done using the polynomial fit to dark energy, 
equation (7). No priors are assumed on w(z). In both panels, the thick solid line shows the best-fit, the light grey contour represents 
the Icr confidence level, and the dark grey contour represents the 2<r confidence level around the best-fit. The horizontal dashed line 
represents ACDM, and the dotted line represents matter density. 




Figure 12. The evolution of w(z) with redshift for fio m = 0.3, using (a) HZT data, and (b) SCP data. The reconstruction is done using 
the polynomial fit to dark energy, equation (7). No priors are assumed on w(z). In both panels, the thick solid line shows the best-fit, the 
light grey contour represents the lcr confidence level, and the dark grey contour represents the la confidence level around the best-fit. 
The horizontal dashed line represents ACDM. 



in this sample is z = 0.86, the behaviour of w beyond this 
redshift cannot be predicted, therefore the apparent flatten- 
ing out of the curve beyond a redshift of one cannot be seen 
in this case. For both these subsets of data, we may repeat 
the exercise using the prior w(z) > —1 (z ^ 0). The results 
obtained for the equation of state, as seen in figures 13(a), 
(b) , are once again commensurate with the results obtained 
earlier for the full sample (figure 6(b)). We may therefore 
conclude from this exercise that subsampling the data does 
not significantly affect our results, and the steep evolution 
of the equation of state of dark energy is not a construct of 
the uneven sampling of the supernovae, but rather, reflects 
the actual nature of dark energy. 



Testing our Ansatz against fiducial dark energy 
models : 

The crucial question of course is whether the recon- 
structed equation of state of dark energy depends upon the 
ansatz which is used in the exercise, i.e. , whether the be- 
haviour of the equation of state merely reflects a bias in the 
ansatz itself. In this section we show how the ansatz per- 
forms in recovering dark energy models whose equation of 
state is known, from simulated data. This ansatz was demon- 
strated to work extremely well when simulations of SNAP 
data were used (Alam et al. 2003). However, simulation of 
SNAP-like data is an optimistic exercise, since data of this 
quality is unlikely to be available in the near future. We 
now demonstrate the accuracy with which the ansatz can 
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Figure 13. The evolution of w(z) with redshift for Qom = 0.3, with the prior w(z) — 1, z 0, using (a) HZT data, and (b) SCP 
data. The reconstruction is done using the polynomial fit to dark energy, equation (7). In both panels, the thick solid line shows the 
best-fit, the light grey contour represents the la confidence level, and the dark grey contour represents the 2a confidence level around 
the best-fit. The hatched region is forbidden by the prior w(z) ^ — 1. The horizontal dashed line represents ACDM. 




Figure 14. Reconstructed equation of state, w(z), for simulated data corresponding to three fiducial dark energy models: (a) quicssence 
with ui = —0.5, (b) generalised Chaplygin gas: p = A/p a , with a = 0.5, wo = —0.9, and (c) w = wo + w\z with too = — 1.2, iui = 1.0. 
f2rjm = 0.3 is assumed and the reconstruction is done using the polynomial fit to dark energy, equation (7). In each panel, the thick solid 
line is the best-fit, the dashed line represents the exact model value, and the light grey contour represents the la confidence level around 
the best fit. 



recover the fiducial background cosmological model if data 
is simulated using present-day observational standards. In 
figures 14 (a), (b), (c), we show how well the ansatz recov- 
ers the equation of state for three fiducial models (assuming 
fiom = 0.3): 

(a) a quiessence dark energy model with a constant 
equation of state: w = —0.5, 

(b) a generalised Chaplygin gas model with p = A/p a : 
with a = 0.5 and the present-day equation of state wo = 
—0.9, which would give rise to an effective equation of state 

W( - Z) = ~\W \ + (1-\1 \)(1+Z)W + ") ' (21) 

and 



(c) a model with a linearly evolving equation of state: 
w(z) = wo + wiz, with w = —1-2, wi = 1. 

(For DE models with w — —1, —2/3, —1/3 the ansatz is 
exact therefore we don't show the results for these cases.) 

We find that in all three cases, the fiducial model lies 
within the 68% confidence limits around the best-fit w(z). 
Based on this result, we claim that within the la error bars, 
the reconstructed equation of state represents the true prop- 
erties of dark energy when we use real data. 

Using other Ansatz : 

It is also important to check whether the results of our 
reconstruction can be replicated using other ansatz such as 
fits to the luminosity distance or the equation of state. Many 
different fits have been suggested in the literature (see for ex- 
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Figure 15. Evolution of w(z) with redshift for Qom = 0.3 with (a) no priors on w(z), and (b) the prior w(z) —1 (z 0). The 
reconstruction is done using Lindcr's fit to the equation of state, equation 22. In both panels, the thick solid line represents the best-fit, 
the dashed line represents ACDM, and the light grey contour represents the Icr confidence level. In the right hand panel, the hatched 
region is forbidden by the prior w(z) ^ — 1. 




Figure 16. Evolution of w(z) with redshift for Qom = 0.3 with (a) no priors on w(z), and (b) the prior w(z) —1 (z 0). The 
reconstruction is done using a four-parameter fit to the equation of state, equation 24. In both panels, the thick solid line represents the 
best-fit, the dashed line represents ACDM, and the light grey contour represents the la confidence level. In the right hand panel, the 
hatched region is forbidden by the prior w(z) ^ — 1. 



ample Huterer & Turner (1999), Saini et al. (2000), Weller 
& Albrecht (2002), Gerke & Efstathiou (2002)). Here we 
choose the fit suggested in Linder (2003) in which the equa- 
tion of state of dark energy is expanded as 



W(z) = + 



WlZ 

l + z 



(22) 



The luminosity distance can therefore be expressed as 



d L {z) 
l + z 



c 

#0 



l + z 



da 



(23) 



where fix = (1 - fi m)a: 3(1+, "° + ''' 1 >exp[3u>i(§ - 1)]. 

We find that for this fit, the errors in the equation of 
state get larger with redshift, however this fit too demon- 
strates that the equation of state of dark energy increases 



rapidly with redshift (figure 15(a)) when no priors are as- 
sumed on the equation of state (EOS). The \ 2 P er degree 
of freedom at the best-fit is Xdof = 1-0298. When the prior 
to(z) > — 1 (z > 0) is invoked, the best-fit EOS remains 
very close to the ACDM model (figure 15(b)). Therefore, 
from this ansatz, we may make the statement that at low 
redshifts, the equation of state of dark energy shows the 
same signs of rising steeply with redshift if no priors are as- 
sumed on the equation of state, thus supporting our earlier 
results. The large errors in the equation of state at redshifts 
of z > 0.5 however make it difficult to make any definitive 
statements about the behaviour of dark energy at high red- 
shifts. 

A limitation of the fit (22) is that it is unable to describe 
very rapid variations in the equation of state. An ansatz 
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which accommodates this possibility has been suggested in 
Bassett et al. (2002) 



w(z) =Wi + ■ 



Wf — Wi 



(24) 



where Wi is the initial equation of state at high redshifts, Zt 
is a transition redshift at which the equation of state falls 
to w(z t ) — (wi + Wf)/2 and A describes the rate of change 
of w(z). 

The resulting luminosity distance has the form: 



d L (z) 
l + z 



c 



r 



dx 



(25) 



where Six = (1 — f2 0m )exp[3 J* 1 (1 + w{z))dz/{\ + z)]. 

The results for the analysis using this fit to the equation 
of state are shown in fig. 16. We find that when the recon- 
struction is done without any priors on the equation of state 
(figure 16(a)), the best fit is remarkably close to the result 
for ansatz (7) (figure 5(b)). The \ 2 per degree of freedom 
at the minimum is Xdof ~ 1-0175 for this fit. The errors in 
this case are somewhat larger, especially at high redshift. If 
we constrain w(z) ^ — 1, then as before, the evolution of the 
equation of state is much slower (figure 16(b)). So the re- 
construction using this ansatz appears to confirm our earlier 
results. 

The above exercises lead us to conclude that our results 
are neither dependent on the nature of the statistical anal- 
ysis nor on the manner in which the SNe data is sampled. 
It therefore appears that dark energy with a steeply evolv- 
ing equation of state provides a compelling alternative to 
a cosmological constant if data are analysed in a prior-free 
manner and the weak energy condition is not imposed by 
hand. 



2.3 Reconstructing dark energy using a new 
Supernova sample 

As this paper was nearing completion, a new dataset consist- 
ing of 23 type la SNe was released by the HZT team (Barris 
et al. 2003). It is clearly important to check whether or 
not these new data points corroborate the findings reported 
in the previous sections. Accordingly, we use a subset of 
200 type la SNe with A v < 0.5 from the 230 SNe sample of 
Tonry et al. (2003), and 22 SNe with A v 0.5 from the new 
sample to obtain a best-fit for our ansatz with fiom = 0.3. 
We then plot the magnitude deviation of our best-fit uni- 
verse from an empty universe with (fiomj^A) — (0.0,0.0) 
in order to illustrate how well our model fits the data (fig- 
ure 17). For clarity, we plot the median values of the data 
points. We obtain medians in redshift bins by requiring that 
each bin has a width of at least 0.25 in logz and contain 
at least 20 SNe. For comparison, we also plot an ACDM 
(fi m, J2a) = (0.3,0.7) model, as well as OCDM and SCDM 
models. From figure 17 we see that our dark energy recon- 
struction is a much better fit to SNe beyond z ~ 0.8 than 
ACDM. At low redshifts (z ~ 0.1) the agreement between 
data and the two models is rather marginal. We now add 
22 of the new supernovae (rejecting one with Ay > 0.5) 
to our existing dataset of 172 supernovae and perform DE 
reconstruction on this new dataset of 194 SNe, assuming 




Figure 17. Literature supernovae (diamonds) shown along with 
median values binned by redshift (large squares). Individual 
points arc shown without error bars for clarity. The solid hor- 
izontal line represents the empty universe with (fJomj^A) = 
(0.0,0.0). The thick solid line represents the magnitude devia- 
tion of our best-fit universe for this data-set from the empty uni- 
verse. The thick dashed line represents ACDM with (Qom, ^a) = 
(0.3,0.7), the dot-dashed and dotted lines represent cosmologies 
with (n 0m ,n A ) = (0.3,0.0) and (1.0,0.0) respectively. 




Figure 18. The evolution of w(z) with redshift for f2o m = 0.3. 
The reconstruction is done using the polynomial fit to dark en- 
ergy, equation (7), for the latest HZT sample of 194 type la SNe. 
The thick solid line shows the best-fit, the light grey contour 
represents the lcr confidence level, and the dark grey contour rep- 
resents the 2cr confidence level around the best-fit. The dashed 
line represents ACDM. No priors are assumed on w(z). 



fiom = 0.3 and no other priors. The resultant figure 18 is 
similar to the figure 5(b), with slightly smaller errors and 
has a best-fit x<Lf = 1-015. The above exercises point to 
the robustness of results reported in previous sections, and 
indicate that evolving dark energy agrees well with the full 
data set containing 194 type la SNe. 
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3 THE ACCELERATING UNIVERSE AND 
THE ENERGY CONDITIONS 

The energy conditions : 

• Strong energy condition: p + 3p > (SEC), 

• Weak energy condition: p^O, p + p ^ (WEC) 

play a vitally important role in our understanding of the 
accelerating universe, both in the context of inflation and 
dark energy. We therefore consider it worthwhile to review 
certain key developments which deepened our understanding 
of these issues. 

In an expanding FRW universe the SEC implies that 
the universe decelerates while the WEC forbids the pressure 
from becoming too negative. Additionally, in the 1960's and 
early 1970's it was noted that energy conditions play a cru- 
cial role in the formulation of the singularity theorems in 
cosmology. Indeed, one of the necessary conditions for the 
existence of an initial/final singularity in big bang cosmology 
is that matter satisfies both the SEC and WEC (Hawking 
& Ellis 1973). 

By the late 1970's it became clear that not all forms of 
matter satisfy the energy conditions. Perhaps the best ex- 
ample of a form of matter which satisfied the weak energy 
condition but violated the strong one is the cosmological 
constant, introduced into cosmology by Einstein in 1917. In 
addition, the vacuum expectation value of the energy mo- 
mentum tensor, (Tifc)vac, which describes quantum effects 
(particle production and vacuum polarization) in an expand- 
ing universe could, in certain cases, violate both WEC and 
SEC (Birrell & Davies 1982; Crib, Mamaev & Mostepanenko 
1980). (For certain space-times, such as de Sitter space, the 
vacuum energy momentum tensor generates a cosmological 
constant since (Tifc) vac = Kgik.) Thus by the late 1970's it 
was well known that neither of the energy conditions could 
be held as being sacrosanct. 5 

The 1980's, as we all know, led to great advances in 
the development of the inflationary paradigm. The inflaton 
field mimics the behaviour of a cosmological constant over 
sufficiently small intervals of time and therefore violates the 
SEC. Early dark energy models were based on inflaton-type 
scalars which coupled minimally to gravity (quintessence). 
Quintessence violates the SEC but respects the WEC. Pre- 
cisely because of the latter property, not any experimentally 
obtained cIl(z) is compatible with quintessence, as empha- 
sized in Sahni & Starobinsky (2000). (The same observation 
holds for H(z), since the latter can be derived from <Il(z) us- 
ing equation (3).) Clearly if observations do indicate that the 
WEC is violated by DE then more general (WEC-violating) 
models for DE should be seriously considered. One example 
of WEC-violating DE is provided by scalar-tensor gravity. 
Scalar-tensor models contain at least two functions of the 
scalar field (dilaton) describing dark energy. As shown in 
Boisseau et al. (2000), these two functions, namely, the 
scalar field potential and its coupling to the Ricci scalar R, 
are sufficiently general to explain any H(z) obtained from 
observations. 

5 The importance of quantum effects to dark energy model build- 
ing has been emphasised in Sahni & Habib (1998); Parker & Raval 
(1999). 



The WEC can also be effectively violated in DE models 
constructed in braneworld cosmology. It has recently been 
shown that such models, with u> e ff < — 1, are in excellent 
agreement with supernova data (Alam & Sahni 2002) . Since 
the field equations in these models are derived from a higher 
dimensional Lagrangian the unusually rapid acceleration of 
the four dimensional universe arises because of the full five 
dimensional theory and not because of matter which con- 
tinues to satisfy the energy conditions and whose density 
remained finite and well behaved at all times (Sahni & 
Shtanov 2003). This behaviour is in contrast to phantom, 
which assumes a conventional 'perfect fluid' form for the 
energy-momentum tensor and therefore contains pathologi- 
cal features such as an energy density which diverges in the 
future and a sound speed which is faster than that of light 
(Caldwell 2002; Sahni & Shtanov 2003). 

The fact that the observed luminosity distance (derived 
from supernova observations) is better fit by dark energy 
violating the WEC than either quintessence or a cosmolog- 
ical constant was first noticed by Caldwell (2002) . Caldwell 
called this 'phantom energy' and showed that larger values of 
Dom (fiom > 0.2) implied increasingly more negative values 
for the equation of state (w < — 1) of phantom. 6 Caldwell's 
results have since been confirmed by larger and better qual- 
ity SNe data sets - for instance Knop et al. (2003) find that, 
in the absence of priors being placed on fiom, the DE equa- 
tion of state has a 99% confidence limit of being < — 1 ! Both 
Caldwell (2002) and Knop et al. (2003) however work under 
the assumption that the equation of state of dark energy is 
unevolving, so that w — constant. 

In this paper we have shown that, suspending the WEC 
prior and allowing the dark energy equation of state to 
evolve brings out dramatically new properties of dark en- 
ergy. Thus the dark energy model which best fits the SNe 
observations has an equation of state which rapidly evolves 
from w(z) < — 1 at present (z — 0) to w(z) ~ at z ~ 1. 
Dark energy therefore appears to have properties which in- 
terpolate between those of dark matter (dust) at early times 
and those of a 'phantom' (w < — 1) at late times. 



4 CONCLUSIONS AND DISCUSSION 

This paper reports the model independent reconstruction of 
the cosmic equation of state of dark energy in which no pri- 
ors are imposed on w(z). In the literature the imposition of 
various priors frequently precedes the analysis of observa- 
tional data sets. Such a procedure is well founded and en- 
tirely justified when priors are dictated by complementary 
information such as orthogonal observations coming from 
different data sets. However, on occasion the use of priors is 
justified on 'theoretical grounds' and in this case one must 



6 Note, however, that we will not consider the theoretical model 
of phantom matter based on a ghost scalar field proposed in this 
paper since, as is well known, it is unstable with respect to parti- 
cle creation (particle + antiparticle of the ghost scalar field plus 
particle + antiparticle of all usual matter fields) and to the loss 
of spatial homogeneity at both quantum and non-linear classical 
levels. 
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be careful so as not to prejudge nature. (Compelling theo- 
retical reasons might well reflect our own particular condi- 
tioning or set of prejudices !) In the case of the analysis of 
type la supernova data, the priors most frequently used have 
been w = constant and w ^ — 1. Both confine DE to within 
a narrow class of models. Moreover, as shown in Maor et al. 
(2002), the imposition of such priors on the cosmic equation 
of state can, on occasion, lead to gross misrepresentations 
of reality. 

In this paper we do not impose any priors on w(z) and 
reconstruct the equation of state of dark energy in a model 
independent manner. In this case our best fit w(z) evolves 
from w < — 1 at z = 0, to w ~ at 0.8 <z< 1.75 (the upper 
limit is set by observations) . This result is robust to changes 
in the value of f2om and remains in place within the broad 
interval 0.1 ^ fiom ^0.5. Our reconstruction clearly favours 
a model of DE whose equation of state metamorphoses from 
w — in the past to w ~ — 1 today. An excellent exam- 
ple of a model which has this property is the Chaplygin 
gas (Kamenshchik, Moschella & Pasquier 2001). However, 
in this model dark energy does not violate the weak energy 
condition (if it was not already violated initially). Our re- 
sults also lend support to the dark energy models discussed 
in Bassett et al. (2002); Corasaniti et al. (2003) in which 
the DE equation of state shows a late-time phase transition. 
An interesting example of an evolving DE model in which 
wo < —1 at present whereas w(z) > — 1 at earlier times is 
provided by the braneworld models (called BRANE1) exam- 
ined in Sahni & Shtanov (2003) which have been shown to 
agree very well with current supernova observations (Alam 
& Sahni 2002). 

It is also conceivable that the observed rapid growth in 
the EOS might characterise 'unified' models of dark mat- 
ter (DM) and dark energy (DE). We end this paper with 
a small speculation on this last possibility. Since the na- 
ture of both DM and DE is currently unknown, it may be 
that a mechanism exists which converts DM (with w = 0) 
into DE (with w ~ —1) in regions with sufficiently high 
density contrast 5p/p ^> 1. (This would happen if, for in- 
stance, the rate of conversion of DM into DE depended upon 
(8p/p) x ,x ^> 1, etc.) Since the conversion of DM to DE is 
confined to high peaks of the density field this process will 
not occur uniformly in the entire universe but will be re- 
stricted to regions occupying a small filling fraction (FF) 
(FF <C 1 for regions with 5p/p 3> 1; see for instance Sheth 
et al. (2003) and references therein). This process could 
commence as early as z ~ 10 — 20 when the first peaks in a 
CDM model collapse. Since DE does not cluster and since 
Pde//0dm grows rapidly as the universe expands, DE from 
high density regions (FF C 1) will spread at the speed of 
light, percolating through the entire universe (FF ~ 1) by 
z ~ 1. Since the creation of DE is tagged to the forma- 
tion of structure, this model may not encounter the 'coinci- 
dence problem' which plagues other scenarios of DE includ- 
ing quintessence. (However this model might have problems 
in producing a sufficiently homogeneous and isotropic dis- 
tribution of dark energy on the largest scales.) The concrete 
mathematical framework for a phenomenological model of 
this kind will be worked out in a companion paper. 

In summary, evolving DE models have been shown to 
satisfy SNe observations just as well (if not better) than the 



cosmological constant. Our best fit equation of state, in the 
absence of any priors, evolves from w(z) < — 1 at z = to 
w(z) ~ at z ~ 1. Indeed, figure 17 shows that our best fit 
EOS is better able to account for the relative brightness of 
supernovae at z > 0.8 than ACDM. However, the evolution 
in w(z) is much weaker if the prior w(z) > — 1 (z > 0) is 
imposed. Due to the paucity of SNe data beyond z = 1.2 
(till date, there is only a single data point beyond z = 1.2, 
SN1999bf at z = 1.75) it is not clear whether w(z) ~ is a 
stable asymptotic value for the reconstructed DE equation 
of state at high redshifts. 7 New supernova data at z > 1 
from ongoing as well as planned surveys (SNAP) combined 
with data from other cosmology experiments (CMB, LSS, 
S-Z survey's, lensing, etc.) are bound to provide important 
insights on the nature of dark energy at high redshifts. Our 
results clearly throw open exciting new possibilities for dark 
energy model building. 
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APPENDIX A: PROPAGATION OF ERRORS 



0.5 



1 .5 



We have seen that the error bars on w(z) for the analysis us- 
ing ansatz (7) are non-monotonic with redshift. Low redshift 
behaviour of the equation of state affects the luminosity dis- 
tance at all higher redshifts, while high redshift behaviour 
effects fewer such distances. This leads to an expectation 
that high-z behaviour of the equation of state should be 
poorly constrained as opposed to the low-z behaviour. This 
seems to contradict the behaviour seen in our figures. To 
investigate if this could be explained by our specific method 
of error analysis we describe the Fisher matrix error bars 
below and show that they are almost identical to what we 
obtain in our method. 

In an analysis which uses an ansatz with n parameters 
Pi, the Fisher information matrix is defined to be 



d 2 L 

dpidpj 



(Al) 



where L = — log£, C being the likelihood . For an unbi- 
ased estimator, the errors on the parameters will follow the 
Cramer- Rao inequality : Apt ^ 1/y/Fu. 

Since the likelihood function is approximately Gaussian 
near the maximum likelihood (ML) point, the covariance 
matrix for a maximum likelihood estimator is given by 



d 2 L 



(A2) 



dpidpj 

The Fisher information matrix is therefore simply the ex- 
pectation value of the inverse of the covariance matrix at 
the ML-point. 

Given the covariance matrix, the error on any cosmo- 
logical quantity Q(pi) is given by : 



2 



i — l j=i+l 



Thus the nature of the errors on a quantity will depend 
essentially on the manner in which it is related to the pa- 
rameters of the system. 

We now consider how errors propagate for different cos- 
mological quantities for the polynomial fit to dark energy 
which we have used for most of the results in this paper : 



H 2 /Hq = QomX 3 



- Aq + A\x + A2X , x = 1 + z 



(A4) 



where Aq = 1 — Qom — Ai — A^. If Oo m is held constant then 
the parameters of the system are {A\, A?). 

We obtain the covariance matrix in {A\,A%) from the 



Figure Al. The deviation of H 2 /Hq from corresponding ACDM 
values over redshift for the ansatz (7). The thick solid line shows 
the best-fit and the light grey contour represents the la confi- 
dence level around the best-fit. The dashed horizontal line denotes 
ACDM. f2om = 0.3 is assumed. 




Figure A2. The variation of the equation of state of dark en- 
ergy w(z) over redshift for the ansatz (7). The thick solid line 
shows the best-fit and the light grey contour represents the la 
confidence level around the best-fit. The dashed horizontal line 
denotes ACDM. Hom = 0.3 is assumed. 



ML analysis, and then using equation (A3), calculate the 
errors on cosmological quantities of interest. For example, 
the errors on the quantity AH 2 — (H 2 — H\ CDM )/H 2 are 
given by : 

<T 2 AH 2(X) = {X- lflCu + 2{X + l)Cl2 + (x + 1) 2 C 22 ] .(A5) 

Although the term C12 is negative we find that a\ H i (x) 
still increases with redshift. This is shown in the figure Al. 
The errors shown are approximately similar to those ob- 
tained in figure 3. 

The corresponding errors on the equation of state can 
be calculated using equations (18) and (A3), and has the 
somewhat more complicated expression : 



\,(x) 



:r 2 [/ 1 2 C 11 +2/ 1 / 2 C 12 + /fC 22 ] 
9[1 - fiom + Ai{x - 1) + A 2 {x 2 - l)] 4 



(A6) 
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where 

h = 1-Oom-A 2 (a;-1) 2 , 

h = 2:r(l-f7 ()m )+AiO-l) 2 , 

and A\, A2 are the mean values of the parameters. Although 
in this case it is difficult to predict the behaviour of error 
bars, after substituting the numerical values we obtain the 
error bars that are shown in figure A2. This figure can be 
compared to the figure 5(b), having almost identical errors. 

This shows that the nature of our error bars is not an 
artifact of our specific method of error analysis. However, 
as shown in figure 15, a two parameter expansion in w(z) 
shows monotonically deteriorating errors in w(z) with the 
redshift, while the expansion in #2 (2) shows errors that im- 
prove with redshift (figure 5(b)). This indicates that the 
nature of error bars might be affected by which quantity is 
being approximated. In the limit of infinite terms in the ex- 
pansion of various quantities all the methods should produce 
identical result. The practical need for truncating these ex- 
pansions make these approximations slightly different from 
each other. More specifically, we require setting of priors 

nz) = ( A7 ) 

i=0 

a n = 0; (n > N v ) (A8) 

where f(z) could be H(z), w(z) or any other physical quan- 
tity and N p is the chosen number of parameters. The non- 
linear priors in the above equation make different finite ex- 
pansions inequivalent. Since we do not know for certain if the 
underlying model for the accelerating expansion involves an 
energy component with negative pressure in a FRW setting 
we are forced to choose one of the alternatives for approxi- 
mations. We hope that with increasingly high quality data 
the effect of such truncations will eventually disappear. 



We find that in both cases, the ansatz recovers the mea- 
sured quantity to within 0.5% accuracy in the redshift range 
important for SNe observations. Thus we find that even for 
fits for which the ansatz does not return exact values, it can 
recover cosmological quantities to a high degree of accuracy. 



APPENDIX B: RECONSTRUCTION OF 
OTHER DARK ENERGY MODELS 

We have seen in the figure 1 that the ansatz (7) works well for 
several physically motivated models of quintessence, Chap- 
lygin gas and SUGRA. In this section we take this exercise 
further and see how well it can reconstruct some of the other 
fits to dark energy known in literature. In figures Bl (a) and 
(b), we show results for simulations using fiom = 0.3 and two 
different fits to the equation of state of dark energy : 

(a) The fit suggested in Linder (2003) : w(z) = wo + 
wi z/(l + z). For this we consider three sets of values in order 
of increasing evolution of w(z) : (a) rao = — 0.8, wi = 0.5, 
(b) wo = — 1.0, wi = 1.0, and (c) too = — 1.2, wi = 2.0, and, 

(b) The non-perturbative w(z) suggested in Corasan- 
iti & Copeland (2003) and Corasaniti et al. (2003), which 
has the parameters Wq (the dark energy equation of state 
today), Wq (the dark energy equation of state at the mat- 
ter dominated epoch), z c (the redshift where equation of 
state changes from Wq to Wq), and A (the width of tran- 
sition). For the simulation we again use three sets of val- 
ues in order of increasing growth rate of w(z) : (a) Wq = 
-0.8, w% = -0.5, z? = 1.5, A = 0.1, (b) w% = -1.0, w% = 
-0.2, z™ = 0.6, A = 0.07, and (c) w Q = -1.2, w% = 
0.1, z™ = 0.11, A = 0.03. 
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(a) (b) 

Figure Bl. The deviation Alog(d^_ffo) /log(d_L-f/o) between actual value and that calculated using the ansatz (7) over redshift for different 
values of parameters for (a) the Linder fit, and (b) the Corasaniti fit for equation of state of dark energy, with f2o m = 0-3. In panel (a), 
the solid line shows the deviation for the Linder fit with wq = —0.8, w\ = 0.5, the dotted line for wq = —1.0, wi = 1.0, and the dot-dashed 
line for wo = —1.2, wi = 2.0. In panel (b), the solid line represents the Corasaniti fit with Wq = —0.8, Wq = —0.5, z™ = 1.5, A = 0.1, the 
dotted line shows w° Q = -1,«>q = -0.2, z™ = 0.6, A = 0.07, and the dot-dashed line, w° Q = -1.2, w™ = 0.1, z™ = 0.11, A = 0.03. The 
dashed horizontal line in both panels represents zero deviation from model values, which is true for ACDM, and ui = — 1/3, w = —2/3 
quicssence models. 



